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Abstract. This paper studies the propagation of intense laser in plasmas in weak 
relativistic region(0.1 < Oq < 0.5) using the quasi-periodic approximation. Effective 
Lorentz factor and density wave effects are calculated in detail. The relativistic 
correction on stimulated Brillouin instability is investigated in the rest parts. The 
coupled dispersion relations of Stimulated Brillouin Scattering(SBS) are obtained and 
investigated numerically. 

PACS numbers: 52.38.-r 
1. Introduction 

The propagation of intense laser through plasmas is an important concern in the laser- 
driven inertial confinement fusion, the laser-plasma accelerators, the x-ray laser and 
other physical problems In General, the amplitude of laser can be described 

by the normalized vector potential {a, 0}, where a = -^s-,<\> = While a C 1, 

the propagating equation of laser is linear. But many non-linear effects will appear 
when a increases , such as density wave effects, relativistic mass correction, etc. While 
a > 1, the problem becomes highly non-linear because of the relativistic effects and it 
can hardly be solved. 

In full ionized plasmas, there are two dominating nonlinear effects on the laser 
propagation: one is the relativistic mass increase of the electron. When the amplitude 
of vector potential a increases, the vibration velocity of electron increases and changes 
the mass of electron to 7m, where the Lorentz factor 7 = , 1 is a coefficient 

determined by the velocity, so non-linear effects appear. 

The other important non-linear effect is the laser's diffraction by electron density 
waves. When the intense laser propagates, electrons in the plasma first vibrate parallel 
to the electric field and perpendicular to the propagating direction of the laser. Then 
the vibrating elections are pushed by the Lorentz force produced by the magnet field 
of the laser. Electron density oscillations parallel to the propagating direction of the 
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laser form density waves. Finally, laser is diffracted by the density waves. As is shown 
in the following, the density fluctuation is proportion to a ■ a, its frequency is twice as 
the laser frequency and its order is the same as relativistic effect's order. So the density 
wave must be considered in the relativistic region. 

In the relativistic region, both non-linear effects cause corrections of the refract 
index of the plasma by modifying the dispersion relation of the laser. The correction is 
very important for the shaping and self-focusing of the intense laser pulse. 

There are many parametric instabilities in plasmas such as Raman instability, 
Brillouin instability and 2uj decay, etc. All these processes will be affected by 
the relativistic mass increase and the density wave effects. Stimulated Brillouin 
Scattering(SBS) is a very important phenomenon because it transfers almost all energy 
to the scattered light and reaches a maximum on backward scattering. In the case of 
backward SBS, it causes laser's reflection and the energy loss. During the implosion of 
ICF, the Brillouin reflectivity can vary from 10% (for short wave incident laser) to more 
than 40% (for long wave laser) In addition, the Brillouin Instability relates closely 
with the filamentation instability. 

In recent papers|H] relativistic mass increase effects are investigated by 
introducing a Lorentz factor. However, because the Lorentz factor varies with the 
electron's jitter velocity, it takes on different values for different non-linear effects and 
should be calculated separately in different cases. Besides, the SBS takes place in all the 
under-dense region of plasmas, while Stimulated Raman Scatter(SRS) can take place 
only in the low density region. Thus the density wave effects should be considered. 

Some authors jTU] [H] already considered the density wave effects with the laser pulse 
propagation without handling the SBS instability. Besides, the density wave equation 
should be altered in the relativistic region (see sec. 2). 

The purpose of this paper is to investigate the Brillouin instability of an intense 
laser with a 2 ~ 0.1 — 0.5. We call this intensity region as Weak- Relativistic Region. 
Which indicates that the nonlinear effects can be processed by the series expansion 
with respect to a 2 . For calculating the Brillouin effect ,we shall use the quasi-periodic 
approximation to handle density wave effects to get non-linear corrections of the laser's 
dispersion relation with higher precision in sec. 2. In sec. 3, we shall investigate the 
Brillouin instability in the relativistic region. 

2. The non-linear effects of the intense laser propagating in Plasmas 

Consider the propagation of a linear-polarized laser in a cold plasma in which the 
electron's thermal velocity is much smaller than their jitter speed in the laser field|12j. 
Suppose the laser is polarized in the x direction and propagates in the z direction with 
a < 1. Assuming the plasma be homogeneous in the transverse direction, we have, from 
the conservation of transverse canonical momentum^, 7/3 x = a x and 7 = y/1 + a 2 , 
where j3 = v/c. 

The laser propagation equation in this case has been derived by P. Sprangle and et 
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al [12] as follows: 

d 2 a d 2 & 2 n 
~ IkH 2 + !h 2 ~ = ( ' 

where n is the density of electrons, n is the local average of n, k p = u p / c and uj p 
is the plasma frequency. 

The plasma should satisfy Vlasov equation. In the cold plasma approximation we 
ignore the thermal-pressure to get the continuity equation and the equation of motion 
from the first two moments of the Vlasov equation: 



gW + ft a(7ft) = 91 _ J__a a 2 (3) 

dct 2 dz dz 27 dz 

The above equations cannot be simplified to an ordinary wave equation as in non- 
relativistic case due to the presence of 7 under derivative operations. The non-linear 
effects must be treated separately for high frequency, low frequency and zero-frequency 
cases to get the correct dispersion relations. 

In the quasi-periodic approximation we assume that the characteristic time of 
evolution of the laser pulse shape is much longer than the laser's period, and a can 
be regarded approximately as a periodic function. Let a = a xcos(cG> t — k z), where x 
is the unit vector in the x direction. Now all of the coefficients and driving forces in 
eq. (j2J) and (j3J) are periodic functions, so its solutions should be periodic far from their 
instability region. 

The approximation is similar to the quasi- stationary approximation used by 
P.Sprangle et al but now the phase velocity of the laser wave in the plasma, uj /ko, 
does not need be close to c. So this approximation can be used for dense plasmas. 

Introduce the phase variable £ = uot — k^z and let 

n = n(£) 

0= 0(0 

After simple calculations, ©becomes 

Uq n d . n — n d<f) Id 



KqC 1 n a£ n oz 27 172; 

Let n = n (l + ip) and use the poisson equation, the linearized equation in ip 



becomes 



d 2 (l^) ,2, 9.1 d , , 



K^-^^M 1 )) (5) 



dcH 2 pT dz y 2~?dz 
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Expand the Lorentz factor 7 and preserve terms lower than fourth, we get the 
nonlinear dispersion relation by an iteration method: 



2 q 1 i.2 



Here the relativistic terms appear with effective Lorentz factors (1 + \afy and 
(1 + |oq) respectively. As we said before, the efficient Lorentz factors differ for different 
modes. The factor ip oc fc 2 , is the frequency-mixing term. The ratio V'oAo increases while 
lu approaches u p . When a 2 , ~ 0.3 and Uq = 2u 2 , we have tp ~ 0.04. 

3. The relativistic correction of Stimulated Brillouin Instability 



Consider a two-dimensional homogenous plasma and ignore the electron's thermal 

JL JL\ Y7 2 — d 2 , d 2 
• 9x' dyh V _L — dx 2 "1" <9?/ 2 ' 



motion in the direction of laser polarization. Let Vj_ = {-§z, 4r}, V 2 = ^ + ' 



(3 = {v x /c, Vy/c, v z /c}, then V 2 = £s + V 2 ± . 

Suppose that the incident laser propagates in z direction and is polarized parallel to 
x axis, the electron's temperature in the plasma is T, then kT <C m v 2 . The conservation 
of the transverse canonical momentum leads to 71^/0 ~ a x . Now the Lorentz force has 
the following form v x (V x a) ^V(a ■ a), 

So the equation of motion for electrons becomes 

^( 7 « + (fl .V)( 7 /J)=V#-±V(...)-V^ (7) 

In stead of P e , we shall introduce the normalized electron thermal-pressure p = 
-, where 6 = kT. 



P e _ (n +n)6 
moc 2 m,(,c 2 

In Brillouin Scattering, the plasma acoustic wave is coupled with the laser wave. 
Let us consider low frequency fluctuations. Ignoring the electron's inertia gives 

{V0}' ow = {— V(a 2 )} Zow 6 V(n low ) 

L J L 2 7 v n m n c 2 v 1 

where {y° w indicates preserving low frequency (acoustic) terms only. 

Substituting this equation into ion's equation of motion, ignoring ion's thermal 
pressure(Tj <C T e ) and using the relation Zni = n we have: 

^ _ lap J [W rn?r + - 

oc 2 t 2 M [ m c 2 n J 

When a < 1, \A + a 2 « 1 + |a 2 — |a 4 . Let a = a L + a s , where a L and a s is for the 
incident laser and the scattered light respectively. 
Now use the relation |a s | <C | , we have 



d 1p s 2w2 Zmoc 2 1 



c"V 2 ^ s = :1 ^-V 2 (a i • a s - ^a L 2 a L • a~ ) low (8) 



dt 2 s T M y " s 2 
- and c s = \ / 

no * V ™ 



where ip s = !3 — and c s = J 4j is the velocity of ion-acoustic wave in the plasma. 
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Similarly, we can write the density fluctuation in the form of if; = — = ipL + ip s an d 
set up the modified propagation equation: 




where is the density wave generated by the laser driving solely. 

The coupled equations (JBJ) and Q are the fundamental equations for the Stimulated 
Brillouin Scattering in the weak relativistic region. There are two differences between 
the present equations and the normal Brillouin equations. First, the factor y/1 + a 2 
have replaced the factor ^a 2 . Second, in the propagation equation, the density wave 
term ipL&s and the relativistic term |a|a s cause a new effect, which is the sideband 
mixing. 

Assume that an acoustic density-perturbation is formed in the plasma: n l ° w = 
noipa cos(u;i — kz). It is coupled with the laser field in the plasma and produces the 
sideband scattered light waves like a s = a + cos[{ojQ + uj)t — (k + k)z] +a_ cos[(uo — u;)t — 
(ko — k)z\. In the linear region(ao <C 1), the two scattered light waves will propagate 
independently and will not affect each other. 

Let u± = Uq ± to, k± = k ± k, in the relativistic region we discussed, the laser will 
drive a density wave ipL whose basic frequency is 2u . This term will cause frequency 
mixing with a + cos(u + t — k + z) and result in a term whose frequency is u_. Similarly, it 
will mix with a„ cos(u;_t — k^z) and result in a term with frequency uj + . Namely, two 
sideband light waves will be coupled to the basic mode. The term a| in |a|a s will act 
in the same way. 

Because of the sideband mixing, a + and a_ must be handled together to get 
cooperative dispersion relations of SBS. 

We consider the simplest sideband mixing effects in order to get rid of the 
unnecessary complexity. Namely, we consider the mixing effects caused by ipi oc 
cos(2cu t — 2k z) and a£ only. Assuming the incident laser propagates parallel to the 
z axis and the scattered light propagate in the y-z plane, the incident wave vector is 
{0, 0, ko}, the two scattered sideband wave vectors are {0, k±, k +k} and {0, — k±, k — k}, 
we can write the incident laser in the form of a^ = aoxcos(c<Jot— koz) , write the scattered 
light as a, = x{a + cos[(u;o + u)t — (k + k)z — k±_y] + a_ cos[(c<Jo — oS)t — (k — k)z + k±y}} 
and write the acoustic density fluctuation as i[) s = ipa cos(ut — kz — k±y). The sideband 
terms in the rhs of the propagation equation 

d 2 d 2 „ 3 

are: 

Sideband terms in ^^a s are 

k 2 ipL&s — ► -^k 2 ipo{a- cos(uj + t — k + z — kj_y) + a + cos(c<j_t — k_z + kj_y)} 

Sideband terms in k 2 ip s a L are 
k 2 

k 2 ip s aLL — > -^-tp a O'o{ cos { l ^+t — k + z — k±y) + cos(u;_t — k_z + k±y)} 
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Finally, sideband terms in — l^a^a, are 
3 3 

— -k 2 & 2 L ai s — > — -k 2 al[a + cos(uj + t — k + z — kj_y) + a_ cos(c<j_t — k^z + kj 

3 3 
~o^l a o a + cosi^-t ~ k- z + kj_y) — -k 2 a a^. cos(cu + t — k + z — k±y) 
8 8 

Matching all the sideband terms, we get: 

D + a + = ^kpipoa- + ^k 2 ip a a - \k 2 ala + - ^k 2 afa- 
D^a^ = 7;kpip a + + ~k 2 ip a a - \k 2 ala_ - ^k 2 ala + 

where D± = — k± — k\ — k 2 . 
a + and a_ can be solved as 

f n - + f k2 p a l)+i\ k l^-l k >D 1 1.2,1, n 

J ^ ~ (D + + |fc2 a 2 )(D _ + | fc 2 a 2 ) _ ( I fc |^_| fc 2 a 2 )22 S^ a 
| _ _ (D+ + |fcgag)+(ifc^o-|fcgag) j , 2 

k °- (D + + |fc2 a 2)(D_ + ffc2ag)-(ifc^ -|^ag) 2 2> 7/; « a 



(10) 



The driving density wave equation is 

— - c 2 V 2 ip s = — V 2 a L • a s - -a?a L • a 

Zitlqc 2 1 3 



2^o(a+ + a_)(l — -a )(A; + fcj_) cos(u;i — kz — k±y) (11) 

Substituting the above expression for a± into it, we get the Stimulated Brillouin 
Scattering's dispersion relation: 

c 2 - c\{k 2 + kl) = ^f\k 2 p (k 2 + kl)(l - 3 -alK 

(D + + D_) + kl^ + \al) 
X ( J D + + |^a 2 )( J D_ + |^a 2 )-[|fc2(^ -|a 2 )P 1 } 

This equation can be solved numerically. Let us consider the backward Brillouin 
scatterjT], when k ~ 2/c and k± = 0. Let u = Qn p c,k = Kk p ,a.nd Qo = uj /uj p , k = ko/k p 
,we normalize the equation (|12|) to 

lc\1 2\ ZtTIq 1 2 / 1 3 2 \ 2 
(S2 — « ) = K (1 «n «n 

v ; M 2 8 o; 



x ■ : ^ : 

SAW 



[fi 2 | - k 2 + |(^o + fog)]" - 4[fifio| - - |(^o - fag 

We adopt the typical parameters ^ = lO -3 ,-^ = 1800 ,the intensity of laser is 
set to a 2 , = 0.32 ,the plasma density is set to Uq = 4u 2 and Uq = 2u 2 respectively. 
The results are plotted on figure 1 and figure 2. In the weak relativistic region, the 
ion acoustic waves become the quasi-mode whose frequency is determined by the laser 
intensity. Correspondingly, the relative value of growing rate and the laser frequency will 
far exceed the ratio c s /c and approach to ^ in the region of the parameter we selected. 
It means that the Brillouin instability will have developed during about several tens to 
a few hundreds periods of the pump laser(about 10~ 13 — 10~ 14 seconds ). 
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Figure 1. the growing rate when n = jUq 



In the region we analyzed, the pure growing modes with Re(u) = appear at large 
k. First, The peak value of the growing rate locates near k = 2k p and has a departure 
due to the relativistic effects. Then, the pure growing modes appear at large k when the 
laser intensity and the plasma density increase. In the weak relativistic region, the pure 
growing modes form a plateau region. The pure growing mode at large k will affect on 
some nonlinear effects like turbulence developing. 

Ignore the relativistic terms and the frequency-mixing terms, the dispersion relation 
(fT2Jl returns to the non-relativistic four-wave dispersion relation of Brillouin scatter. 
We can solve the non-relativistic dispersion relation and compare the result with the 
relativistic result. All the results are plotted in figure 3. In which we can find the effects 
of relativity and frequency-mixing are: (l)move the peak value; (2)reduce the growing 
rate of instability. The reduction is larger in the plateau region of pure growing modes. 
So the Brillouin reflection rate will be reduced by the effects. 

4. Conclusion 

The present paper applied the quasi-periodic approximate method for the stable laser 
propagation in weak relativistic plasmas. The Stimulated Brillouin Scatter in relativistic 
region is investigated in detail. The growing rates are re-calculated numerically and the 
sideband-mixing effects are considered. The Brillouin instability is reduced by these 
effects. For the plateau region of pure growing modes, the reduction is more evidence. 
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Figure 3. the comparison of relativistic and non-relativistic growing rates. 
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